Limits of dark energy measurements from CMB lensing-ISW-galaxy count correlations 
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I discuss several issues that arise when trying to constrain the dark energy equation of state 
using correlations of the integrated Sachs- Wolfe effect with galaxy counts and lensing of the cosmic 
microwave background. These techniques are complementary to others such as galaxy shear surveys, 
and can use data that will already be obtained from currently planned observations. In regimes where 
cosmic variance and shot noise are the dominant sources of error, constraints could be made on the 
mean equation of state to ±0.33 and its first derivative to ±1.0. Perhaps more interesting is that 
the determination of dark energy parameters by these types of experiments depends strongly on the 
presence or absence of perturbations in the dark energy fluid. 

PACS numbers: 98.80.Es 



I. INTRODUCTION 

Understanding cosmic acceleration is one of the largest 
problems facing physics today. So far the most direct 
measurements of acceleration have come from distance- 
redshift measurements 0, 0, 0- This acceleration is 
thought to be due to the effect of dark energy, a new 
form of energy density that dominates the current uni- 
verse. However, a universe with dark energy exhibits a 
different evolution for density perturbations, and this has 
observable consequences. 

In this paper, I will discuss how well certain measure- 
ments of growth can constrain the expansion history of 
the universe, and thus dark energy. In particular I will 
discuss the lensing of the cosmic microwave background 
(CMB) by intervening structure combined with measure- 
ments of galaxy number density correlations from sur- 
veys. These types of measurements are attractive since 
they can be obtained essentially "for free" from experi- 
ments already planned. I will also discuss the sensitivity 
of these measurements to the presence of perturbations 
in the dark energy fluid, which can have an important 
effect on the results Q. In this paper "lensing" refers 
specifically to lensing of the CMB by foreground struc- 
ture, rather than the extremely productive research area 
of studying the shape distortion of background galaxies 
due to lensing from foreground galaxies (covered exten- 
sively in |5|,|6|). 



II. BACKGROUND THEORY 



A. The ISW effect and dark energy 



blueshifts and redshifts will cancel each other out, leav- 
ing no net effect. However, if the potential does change 
over time there may be some overall change in each pho- 
ton's wavelength, and hence the observed temperature. 
This is the integrated Sachs- Wolfe (ISW) effect, with the 
change in temperature AT ISW observed in a direction n 
expressed simply as 



AT ISW (n) = -2 / dD$(x(n),D), 



(1) 



where the gravitational potential is written as a function 
of position x and lookback distance D (used as a proxy 
for conformal time). 

For a flat, matter-dominated universe, the potential 
remains constant over time even though the density per- 
turbations themselves do not. This means that any ISW 
effect originates either from early times, when the den- 
sity of radiation was still significant enough to affect the 
expansion rate, or from late times when dark energy be- 
came dominant. Thus if the late ISW signal can be sep- 
arated out from others, it provides a clean measurement 
of dark energy. 

Normally the ISW effect itself is buried by the primary 
temperature anisotropies from the LSS. However, the pri- 
mary anisotropics are set at the last scattering surface, 
at a completely different epoch and at different length 
scales than those of the structure growth responsible for 
the ISW effect. The primary anisotropy should thereby 
be uncorrelated with the ISW effect and other measure- 
ments of growth. Thus cross-correlating other measure- 
ments with CMB temperature maps can be a useful tool 
for bringing out information about growth, as has been 
discussed recently by 0, and demonstrated by . 



As photons travel from the last scattering surface 
(LSS) to us, they fall into and climb out of potential wells 
that lie along their path. If the gravitational potential 
<J> does not change with time, then the accompanying 
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B. Lensing correlations 

The first measurement I will consider correlating with 
the temperature map is a measurement of gravitational 
lensing of the CMB. The microwave background is grav- 
itationally lensed by matter that lies between us and the 
LSS. The map of photon deflection angles over the sky 
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can be written as the gradient of a scalar field <fi called 
the projected potential, which depends on the 3D gravi- 
tational potential $ as 



(n) = -2 



dD D ™ S D ^ $>(x(h),D), (2) 



where -Dlss is the distance to the last scattering surface. 

Just as with CMB temperature maps, the map of the 
projected potential can be decomposed into spherical 
harmonics. Then a two-point function of the modes can 
be put together to construct an angular power spectrum. 
The expression for the power spectrum can be written as 
a line-of-sight integral over the 3D gravitational potential 
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dD D <& 2 (fc, D) 



Ds-D 
DD S 



k=e 



(3) 

While this equation is only exact in the flat-sky (large t) 
limit, the important feature is that it captures the physics 
of how the angular power spectrum depends on a line-of- 
sight integral of the gravitational potential multiplied by 
a kind of window function and the primordial (i.e. the 
growth function has been separated out) power spectrum 
Pg(k). The exact expressions and further details can 
be found in [lj}; they are used for the computations in 
this paper for multipole moments where the difference is 
important. However, for the remainder of this section 
1 will write only the Limber-approximated integrals for 
clarity. The cross correlation between the temperature 
and lensing is due to the presence of the ISW effect in 
the temperature, and thus its power spectrum has the 
simple flat-sky form 
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C. Galaxy correlations 
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FIG. 1: Noise levels for CMB lensing reconstruction. The 
dark curve is the lensing power spectrum for a typical cos- 
mology. The different dotted lines are noise levels for CMB 
temperature experiments with a 4' beam and 0.1 fj,K noise per 
pixel (dotted line) or 0.01 \iK noise per pixel or less (dashed 
line). The lower dashed line is for an E-mode CMB polariza- 
tion measurement. 



and with the lensing potential power spectrum 



J dDD$ 2 (k,D)n g (D) 



DD S 



(7) 

Again, in the calculations described later these Limber- 
approximated integrals are used only for high multipole 
moments. The exact expressions are used for low multi- 
pole moments (I < 100). 



The second type of observation I will consider is count- 
ing galaxies projected on the sky. On large scales, fluc- 
tuations in the number density of galaxies should track 
fluctuations in the gravitational potential (possibly with 
some bias). From the map of number density over the 
sky, we can obtain the auto-correlation power spectrum 



r<99 



dDD^ 2 (k,D)n 2 a {D)Pl{k)\ k __ 



(5) 



which depends on the potential as viewed by some win- 
dow function n g (D) describing the distribution of galax- 
ies that are actually observed. The galaxy spectrum 
should be correlated both with the ISW part of the tem- 
perature spectrum 
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(6) 



III. CALCULATING THE SENSITIVITY TO 
DARK ENERGY 

The power spectra themselves are computed numeri- 
cally using the techniques described in the previous sec- 
tion and a version of the CMBFAST code |ll| which I mod- 
ified to output lensing spectra and other information. 

The basis for all the analysis is the Fisher informa- 
tion matrix. Given a fiducial model, the Fisher matrix 
describes how sensitive the model is to changes in its pa- 
rameters. First, all the power spectra are put together 
into a covariance matrix 



G, 



TT 



c, 



c t = 



J* " Cf + Nf Cf 9 | . IN) 
T 9 Cf 9 C 99 + N 9 



3 



Experimental noise can be modeled as a contribution to 
the diagonal terms of this matrix. In our case I take the 
noise contribution to the primary temperature anisotropy 
Nj to be negligible. However, even a small amount of 
noise is significant when reconstructing lcnsing informa- 
tion, and I take this into account using the derivations 
for lensing error as a function of temperature noise found 
in |T^ and shown in Fig. ^ I take the primary contribu- 
tion to the noise for galaxy surveys Nf to be shot noise 
(bias for each survey is treated as a free model parame- 
ter). Some other possible noise sources are discussed in 
section lill Bl 

When I consider combining three galaxy surveys at 
different redshifts, each galaxy survey gets its own row 
and column, resulting in a 5 x 5 covariance matrix. The 
Fisher matrix is then constructed as a sum over multipolc 
moments, 



F, 



dsi 



ds 



(9) 



where the Si are labels for the actual model parameters. 

The model parameters actually varied are the dark en- 
ergy parameters discussed in the next section, plus an 
angular scale parameter £a, baryon density flbh 2 , mat- 
ter density fl m h 2 , primordial power amplitude As and 
tilt ns, optical depth r, and a bias parameter hi for 
each galaxy survey. All models were constrained to be 
flat. I investigated two fiducial models; both had pa- 
rameters equivalent to a modern concordance cosmol- 
ogy (fi d . e . = 0.75, fl b = 0.04, Q cdm = 0.21, h = 0.65, 
t = 0.1, ns = 0.9). The difference between the two was 
in the dark energy evolution; one model had a "pure A" 
equation of state of w = —1, whereas the other used a 
quintessence-type model with w = —0.9 today running 
smoothly to w = — 1 at high redshift. The latter has a 
somewhat enhanced sensitivity to dark energy parame- 
ters due to the evolution of dark energy at low redshifts, 
and is used for reporting the final constraints on dark en- 
ergy parameters in section llVl Using the "pure A" model 
instead degrades the dark energy parameter constraints 
by about 20% overall. 



A. Choosing a form for dark energy 

Different experiments have differing sensitivities to 
dark energy. With no knowledge of the fundamental 
physics behind cosmic acceleration, there is little rea- 
son to favor one functional form for the equation of state 
w(z) over another. In that light, principal mode analysis 
is useful for revealing what the experiments can actually 
constrain. In such an analysis one chooses some basis set 
of eigenfunctions (cut off to a finite set size), and then 
expresses the fundamental modes for each experiment as 
combinations of basis functions. In our case the basis 
functions will simply be a set of 25 boxcar functions cov- 
ering the redshift range from z = to z = 5. 
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FIG. 2: Three dark energy eigenvectors for an ideal ISW- 
lensing- galaxy count correlation experiment. The left part of 
each plot shows the 25 w(z) parameters, and on the right are 
the 9 other parameters (in order): angular scale I a, baryon 
density Qbh 2 , matter density Q m h 2 , optical depth r, primor- 
dial power spectrum tilt ns, primordial power amplitude As, 
and the bias parameters for each galaxy survey 61, 62, 63- The 
error in each eigenmode for a cosmic-variance limited experi- 
ment is shown for each mode. 

A few example eigenvectors are shown in Fig. [21 In 
general, dark energy information is mixed with other cos- 
mological parameters so that the eigenmodes are compli- 
cated combinations of several parameters. However, with 
the Fisher matrix in hand not only can one perform prin- 
cipal mode analysis, one can also marginalize over some 
parameters and switch to more familiar dark energy pa- 
rameters with simple matrix operations. For example, it 
is simple to switch from the 25 w(z) parameters to the 
"u> a " parametrization 0,0, which takes the form 



w(z) = w + w a 



1 



(10) 



The full eigenmode analysis shows that even for ideal 
experiments, no more than two or three dark energy pa- 
rameters will be well-constrained, so for the rest of this 
paper the "w a " parameterization will be assumed. 



B. Galaxy survey characteristics 

My intent is to approximate the behavior of galaxy 
surveys which do not include spectra, but are able to 
obtain approximate photometric redshifts through color 
information. I therefore consider the situation where one 
has a large number of galaxies that can be assigned to one 



4 



1 


1 1 1 1 lllj 1 1 


1 1 inij i 


i 1 1 inij i 


i 1 1 1 nij i 


1 1 1 in 


- 


(5w a 








i 


- 




* 


* 




- 
















(5w 






















1 


1 i i 


1 i 


1 i 


1 i 






' 10 7 10 8 10 9 10 10 10 11 10 1S 
number of galaxies 

FIG. 3: Error in wo and w a as a function of total galaxy 
number with lensing noise fixed at the "< 10 nK E-mode" 
level from Fig. The solid lines with squares are for a sur- 
vey with three redshift bins, dashed lines with triangles for a 
survey with one high redshift bin, and dotted lines with stars 
for one with a single low redshift bin. 



of three redshift bins, centered around z = 0.5, z = 1.0, 
and z = 1.5. Each of these bins is approximated as a 
rounded boxcar function similar to the technique of Hu 
and Scranton [15j . This approximation is best for survey 
slices that are not severely magnitude or volume limited. 
For those cases the overall effect is to smear out the edges 
of the bins, weakening the advantage of having several 
redshift bins. 

Galaxy surveys may contain numerous systematics. 
Stellar contamination should not be a big problem since 
for faint surveys the majority of observed objects are 
galaxies anyway. More worrisome is that the analysis 
presumes knowledge of the redshift distribution of galax- 
ies in the sample, which of course could be in error. We 
can estimate how small the overall redshift calibration 
error needs to be by converting the redshift error into 
an error in the angular scale of galaxy correlations, and 
examining the change in the correlation power spectrum 
compared to the shot noise. Because of the large an- 
gular scales we're interested in, it turns out that even 
for galaxy surveys with 10 12 galaxies (corresponding to 
roughly 10 4 per square arcminute), shot noise at these 
scales is still large enough that an overall redshift cal- 
ibration error of up to a percent is tolerable. For the 
remainder of this paper I will assume any systematics in 
the galaxy correlations are below the shot noise level, but 
new experiments may well reveal important new system- 
atics not considered here. 



IV. RESULTS 

In the absence of noise, over 40% of the improvement 
in dark energy constraints gained by adding CMB lensing 
data to the temperature information comes from the Cj^ 
cross-correlation channel. Even though in the ideal case 
the noise in the CMB and lensing potential anisotropics 
is uncorrelated, the inversion of the covariance matrix in 
equation[5]means that noise in the auto-correlation spec- 
tra will find its way into the cross-correlation channel. 
This then means that the noise in the reconstructed 
spectrum will be the limiting factor as to how well CMB 
lensing can constrain dark energy. 

In Fig. I have calculated the noise in the recon- 
structed Cf^ lensing power spectrum for several different 
types of experiments. At small enough temperature noise 
levels (a bit above 10 nK per pixel for 4' pixels) all of the 
noise in the lensing Cf^ is in fact coming from cosmic 
variance in the temperature Cj T . Note that even at this 
limit, the noise in the lensing spectrum is still roughly the 
same order of magnitude as the signal. Part of the low 
signal-to-noise ratio appears to be related to the choice of 
a red tilt (ris < 1) in the primordial spectrum of the fidu- 
cial model. This reduces CMB power at the small scales 
from which the lensing signal is reconstructed, leading to 
more relative noise in the lensing power spectrum. The 
total noise in the lensing power spectrum can also be re- 
duced further by combining several polarization modes 
as per 12], which results in an improvement by a factor 
of a few, so it may not be necessary to drive the detec- 
tor temperature all the way down to WnK to get the 
desired precision. In any case cosmic variance is the ul- 
timate limit that I will consider here, and we shall see 
that it is a fairly restrictive one, at least as far as w(z) 
constraints are concerned. 

How about cross-correlating galaxies with the ISW ef- 
fect in order to augment our information about dark 
energy? This can provide interesting results, especially 
with the possibility of adding even very limited redshift 
information into the mix. I show the results in Fig. [3] 
Galaxy information, especially that coming from high 
redshift galaxies, can tighten the dark energy constraints 
by more than a factor of two. 

Finally, the constraint contours in the wo-w a plane 
are shown in Fig.^J These are constraints made on w(z) 
simultaneously with constraints on all the other cosmo- 
logical parameters described above. In the absence of 
any other experiments, CMB lensing and galaxy counts 
combined can realistically constrain (contour b) wq to a 
precision of about ±0.33 when all other parameters are 
marginalized over. 



A. Dark energy perturbations 

The above results were obtained assuming the absence 
of perturbations in the dark energy fluid. A true cosmo- 
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FIG. 4: Error contours (68.3%) in the wo~w a plane. From 
innermost to outermost they represent (a) a "perfect" CMB 
lensing- galaxy count cross correlation measurement where 
cosmic variance is the only limitation, (b) a more realistic 
measurement with galaxy density of ~ 10 2 arcmin" 2 and the 
"< 10 nK E-mode" lensing noise of Fig. Q (c) an ISW-lensing 
only experiment (no galaxy counts). 



logical constant has w = — 1 and is perfectly smooth with 
no perturbations. However, many quintessence models 
are built with some sort of scalar field which in gen- 
eral can have its own density fluctuations. These den- 
sity fluctuations become important at late times when 
the structure formation responsible for the 1SW effect is 
occurring, and thus significantly affect the above results. 
This effect is shown in Fig. |SJ Note that including per- 
turbations mostly increases the uncertainties along the 
degeneracy direction, especially for experiments limited 
only by cosmic variance, with the result that limits on 
both dark energy parameters are degraded by nearly a 
factor of three. 



V. CONCLUSION 

Cross-correlating the ISW effect with CMB lensing and 
galaxy counts can in principle place limits on dark en- 
ergy; ideally (if cosmic variance were the only limita- 
tion) one could measure wq to ±0.093 and w a to ±0.32, 
though in practice realistic limits are probably worse by 
about a factor of three. These measurements depend on 
dark energy's effect on the growth of structure, and thus 
have a dependence on the equation of state quite differ- 
ent from distance-redshift measurements. Observation of 
such growth effects would bolster the case for dark energy 




FIG. 5: The effect of dark energy perturbations on error con- 
tours (68.3%) in the wo~w a plane. Contours (a) and (b) are 
as before in Fig. 0] Contours (c) and (d) show the effect of in- 
cluding dark energy perturbations for experiments considered 
in (a) and (b), respectively. 



as the source of cosmic acceleration. Also, these measure- 
ments are strongly sensitive to perturbations in the dark 
energy fluid (which do not affect distance observations) 
and thus may ultimately be more useful as a measure- 
ment of perturbations than as a precise determination of 
the equation of state. 

There is certainly room for improvement for these 
kinds of measurements. The amount of lensing being 
measured is small, so cosmic variance and the finite res- 
olution of experiments become important. In this work I 
have used the CMB temperature map to reconstruct the 
lensing spectrum; this method essentially makes use of 
four-point correlations in the original temperature map. 
The three-point correlation function (known as the bis- 
pectrum) also contains contributions from lensing which 
are known to be significant and potentially less noisy, 
although how well dark energy information can be ex- 
tracted from realistic data is still under investigation. 

The scales on which galaxy densities correlate with the 
ISW effect are very large (degree scale or larger), so the 
deviations from a smooth background are small and even 
for large numbers of galaxies shot noise is significant. One 
way of overcoming this problem is simply to use more in- 
formation from galaxies than merely their number den- 
sity. Dedicated telescopes that can find galaxy shapes 
and redshifts are predicted to put good constraints on 
dark energy [f|. In the future a wide array of comple- 
mentary observations will be available to determine the 
nature of cosmic acceleration, each with their own sensi- 
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tivity and limitations. 
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